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Q^ ' We study by Monte Carlo techniques the evolution of finite two-dimensional 
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A periodically modulated bistable system in the presence of noise exhibits the 
phenomenon of stochastic resonance (SR) [1-3]. Plotting the correlation a between 
the modulation signal and the response of the system, versus the noise intensity 
we usually obtain a strong peak. This characteristic picture is accounted in the 
literature as the phenomenon of SR. A brief review of the theoretical aspects is 
given in [4]. Experimental evidence for the SR phenomenon was found in analog 
simulations with proper electronic circuits [3,5,6], laser systems operating in multi- 
stable conditions [7], electron paramagnetic resonance [8,9] and in a free standing 
magnetoelastic ribbon [10]. 

In the present paper we report on the possibility of obtaining the SR in a 
bidimensional Ising system. The phenomenon of SR was already studied in globally 
coupled two-state systems [11,12]. Due to the fact that the spins in the Ising model 
can be considered as coupled bistable elements, the proposed system is special case 
of the problem treated in [11]. In contrast with all earlier works, we do not consider 
any external stochastic forces, just the thermal fluctuations in the system. 

Considering an Ising system on a finite square lattice at zero thermodynamic 
temperature, the free-energy curve versus the magnetization will have a double well 
form [13]. An external oscillating magnetic field will modulate the two minimums 
in antiphase. The effect of a positive temperature in this system can be consid- 
ered as a stochastic driving force, and the magnetization in function of the time 
{m{t) = J2i Si) as the response function of the system. In this way all the necessary 
conditions are satisfied for the SR phenomenon. Due to the fact that the noise 
intensity (thermal fluctuations) is temperature dependent, the characteristic peak 
must be observed for a Tj. resonance temperature in the plot of the correlation (a) 
versus the temperature (T). 

It worth mentioning that recently the Ising system in oscillating fields was con- 
sidered by computer simulations, studying the hysteretic response of the system 
[14]. 

We studied the SR in the proposed system by a Monte Carlo method using the 



well-known Metropolis algorithm [15]. The Hamiltonian of the problem is: 

H=-JJ2S^S] + f,BB{t)J2S^, (1) 

where the sum is refering to all nearest-neighbours, S*/ = ±1, /i^ is the Bohr mag- 
neton and B{t) the external magnetic field. We will consider B{t) in a harmonic 

form: 

27r 
B{t) = A.s\n{-t) (2) 

The time-scale was chosen in a convenient form, seting the unit-time interval equal 
to the average characteristic time (r) necessary for the flip of a spin. We have 
taken this r time interval as a constant, and thus independent of the temperature. 
Although this assumption is just a working hypothesis, we expect usefull qualitative 
results. 

Our Monte Carlo (MC) simulations were performed on square lattices with 
N X N spins, considering the value of A^ up to 200. One MC step is defined as 
N X N trials of changing spin orientations and corresponds to a r time-interval. 
(The period (P) of the oscillating magnetic field is also given in this r units.) The 
amplitude (^4) of the magnetic field is considered already multiplied with ^, and 
thus it has the dimensionality of the temperature {k is the Boltzmann constant). 
The temperature is given in arbitrary units. The critical temperature of the infinite 
system 

T, = 2.2692.. .p (3) 

is always considered as 100 units. Starting the system from a completely random 
configuration, to approach the dynamic thermodynamic equillibrum we considered 
5000 MC steps. The correlation function between the driving field and the magnetic 
response of the system 

1 " 
a =< Bit) ■ m{t) >= - • E B{ti) ■ m(t,), (4) 

^ j=i 

was studied during n = 5000 extra iterations. (The averaging in the expression of a 
is in function of time, and ti = t ■ i) The correlation (4) was studied in function of: 



• the temperature (T) 

• lattice size (A^) 

• amplitude of the magnetic field {A) 

• period of the magnetic field (P) 

Our results are summarized in Fig. 1-4. As we expected for a fixed lattice size and 
a given oscillating magnetic field the curve a versus T exhibits the characteristic peak 
of the SR. Considering A = 10, P = 50 and A^ = 5 a generic result is ploted in Fig.l. 
From Fig. 2 we conclude that the location of the peak is not significantly infiuenced 
by the amplitude of the magnetic field. As expected the amplitude infiuenced only 
the shape of the peak, and its height is stongly increasing with the amplitude. (Our 
results would suggest an exponential type variation.) For small lattices (A^ < 10) 
this resonance temperature (T^) exhibits a strong dependence of the lattice size, 
and converge to a T^ limit value for big lattices (A^ > 100): 

Tr = liniN^ooT^ (5) 

In Fig. 3 we present the results for three choices of P. One can immediately observe 
that the T^ resonance temperature is dependent on the period of the magnetic field. 
As we illustrated in Fig. 4 this dependence can be approximated with the curve: 

Tr = T^ + 557.264 ■ p-o-6335 (g) 

In the limit of high frequencies the resonance temperature is tending to infinity 
and in the limit of small ones, Tr is tending to Tc. Due to the fact that the real 
experimental conditions work in the very low frequency limit (the periods in Fig. 4 
are given in units of r!), one would expect this phenomenon detected at T^. 

Due to the sensibility of T^ on the size of the system (Fig. 3 for small A^) an 
interesting subject would be the study of this effect on materials containing small 
magnetic domains. 

The theory of SR has been previously developed in the context of classical sta- 
tistical physics using linear response theory and the fiuctuation-dissipation relations 



[16,17]. In this sense, by calculating the susceptibility of our model in the absence 
of the periodic magnetic field and studying the response of the system with the 
linear response theory one could perform analytical studies of the problem. We 
consider that such studies would be important both for statistical physics and for 
the phenomenon of SR. 

From computational point of view, a much more detailed study on the phe- 
nomenon, and the consideration of the real 3D case would also be of interest. 

In conclusion, in this letter we studied by MC techniques the phenomenon of 
SR in finite Ising systems considered in a periodic magnetic field. The thermal 
fluctuations were considered as a stochastic force, and the resonance was detected 
for a Tr resonance temperature. The dependence of T^ against the characteristics 
of the magnetic field and the lattice size were investigated. 

The MC simulations were performed on the computers of the Theoretical Physics 
Department (SENTEF) from the University of Bergen (Norway). I am greatfull to 
L. Csernai for his continuous help. 
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Figure Captions 



Fig. 1. Characteristic peak for the SR phenomenon in the plot of the correlation 
{a =< B{t) ■ M{t) >) versus the systems temperature [Tc = 100, A = 10, P = 50 
and N = 5). 

Fig. 2. Shape of the peak for three different values of the magnetic fields 
amplitude {A = 5, 10 and 20) (T^ = 100, P = 50 and TV = 5). 



Fig. 3. Dependence of the resonance temperature versus the lattice linear size 
for three different values of the period (20, 50 and 200) {Tc = 100 and A = 10). 



Fig. 4. Variation of the resonance temperature (T^ — Tc) for big lattices (200 x 
200) against the period of the magnetic field. The best-fit curve indicate: Tr — Tc = 
557.264 ■ p-o-6335^ 



